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SECTION I 


INTRODUCTION 


Composite laminates are increasingly being used in the design of 
structural components of flight vehicles. To evolve an efficient design 
of components, it is necessary to know behavior of the laminates under 
various combinations of loads. This behavior can be determined experi- 
mentally but this is both expensive and inconvenient. The alternative 
to the purely experimental technique is to determine the analytical 


response of laminates based upon the mechanical properties of unidirec- 


tional laminas and to conduct some limited experimental verification. 


Analytical technique relating the properties of the unidirectional 
laminas to those of the laminates is documented in References 1 and 2. 

The concepts developed in References 1 and 2 form the basis of a 
NOnlinear LAminate STrength (NOLAST) analysis computer program which 
is the subject of this report. The NOLAST program predicts the rsponse 
of multidirectional laminates subjected to in-plane loads. 

A general description of the program is included in Section II 
while Sections III and IV contain user's instructions and program 
modifications along with sample problems. The Fortran listing of 
the program and examples of output appear in Appendices A and B 


respectively. 


SECTION II 


PROGRAM DESCRIPTION 


The program NOLAST is designed to analyze multidirectional laminates 
subjected to monotonically increasing proportional in-plane loads to 
failure. In the prediction process, it uses the following: 


1. The incremental constitutive relationship for a unidirectional 


lamina 


de (c,) do, (i, j = 1, 2, 6) (1) 


¢” Say j 
where dey» dos Sij (e,) are stress and strain components and elements of 
the compliance matrix. For accurate determination of S44 etc.. experi- 
mentally obtained basic stress-strain data (uniaxial tension and compres- 
sion along awa transverse to the material axes and longitudinal shear) of 
the unidirectional lamina are represented by piecewise cubic spline inter- 
polation functions to yield smooth composite stress-strain curves. The 
spacing of the stress-strain points of the experimental data can be chosen 
arbitrarily. 

2. In paragraph 1, elements of the compliance matrix, Say (ey), are 
functions of strains (e,) resulting from simple load conditions. In case 
the stress-strain curves of the lamina are highly nonlinear, the use of 
strain components under biaxial stresses would have significant effect on 
the predicted behavior of the laminate. The presence of a compressive 
stress in the transverse dire tion in combination with a tensile stress in 


the longitudinal direction is likely to result in the reduction of the 
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apparent longitudinal tensile tangenet modulus. The reverse situation 
would occur if the transverse stress were tensile. The quantitative 
influence on the tangent modulus relative to one material axis due to 
the presence of stresses in the other material axis is not know. This 
information is not available due to the lack of experimental data under 
biaxial stress states. Consequently the use of components of strains 
under biaxial stress fields to determine S5j (e,) is erroneous. For 
the determination of the tangent moduli in this program the effect of 
the presence of transverse stresses (or longitudinal) is allowed for. 
For this purpose it is assumed that simple equivalent strain increments 


can be computed from the following Equations. 


* 


ac,| 


= de, /(1-v, 5B) (2) 
Eq. 


= de,/(1-v,,/B) (3) 
Eq. 2 21 


where 


Yy2 ™ Major Poisson's Ratio and 


B= do,/do, 


3 In order to compute laminate strain increments, Equation 24 of 


Reference 1 is used. This Equation is 


[dee] = (a}7? fan] (4) 


where [de°] = Laminate Strain Increments 

tay = Laminate Compliance Matrix 

[dN] = Increments of Stress Resultants 
The compliance matrix ray in Equation 4 represents the average tompliance 
properties during the application of (nt+1)-th load increment. However, 
these properties are not known beforehand. This difficulty is overcome by 
using the elastic properties corresponding to those at the end of the n-th 
load increment to compuce laminate strain increments. These are then used 
to compute current stresses and strains in the plies. Based upon the cur- 
rent strains in the plies, average elastic properties are determined. A 
new compliance matrix is computed and laminate strain increments determined. 
This procedure is repeated until the difference between two laminate strain 
increments in two consecutive cycles is smaller than the prescribed error 
bound. This predictor-corrector and iterative procedure renders the method 
of analysis practically independent of the size of the load increment. 

4. Incremental loading of the laminate cannot continue indefinitely 
without affecting the load carrying capability of the plies. To determine 
the onset of degradation of the plies, a criterion is used. The criterion 
used in this program (Equations 21, and 22 of Reference 2) is a function of 
both stress and strain states. Specialized for plane stress conditions, it 


can be written as 


a | 
Ky if o,de,] “= 1 (i = 1, 2, 6) (5) 
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where 


k= cf o,dey] * (i= 1, 2, 6) (6) 


In Equations 5 and 65 He te are current and ultimate normal (tensile or 
compressive) and shear components of strains. 
In plywise application of Equation 5, the program NOLAST allows for 


two modes of failure i.e. fiber failure and matrix failure modes. To 


ee 


distinguish between the two modes, it is assumed that if 


K [ O,de ] K f o,de Zz l 7 


the failure mode of the ply corresponds to that of the fiber. When the 
ply reaches a degradation state, the ply is assumed to unload while the 
laminate loads are maintained. In the case of a matrix failure mode, 
both transverse and shear loads of the affected lumina are set to zero. 
If the failure corresponds to that of the fiber, total unloading of the 
affected plies is assumed. 

5. Strain components corresponding to biaxial stress states are 
used to satisfy Equations 5 to 7, while equivalent strain components are 
used to determine tangent moduli. In the absence of experimental data, 
it is assumed that these strain components cannot be in excess of the 
ultimate strains obtained under simple load conditions. For this 


reason, these two type of strain components are checked at the end of 


application of each load increme-t. If the difference between ultimate 


experimental strain components and equivalent and/or biaxial strain com- 
ponents is less than a prescribed value, the ply is assumed to have failed. 
The mode of failure would depend upon the nature of the strain component. 
The exceedance of longitudinal strains results in the fiber failure mode 
while the others cause a matrix mode of failure. 

6. Once the ply failure and its associated mode have been estab- 
lished (Point B in Figure 1), there still remains a question of the 
influence of the failed ply upon the laminate. Some of the possible 
responses of the laminate subsequent to the initial ply failure are shown 
in Figure 1. They are the following: : 

a. The affected ply unloads completely at B giving rise to a : 
instantaneous increase in strain BB, or a drop in stress BB, depending | 
upon the nature loading process. Under constant loading rate, the path | 
Bi, would be obtained, while the path BB, would result under constant | 
displacement application. After By or Bo» the stress strain response ; 
would be along BAA, or BAA: : 

b. An other extreme response corresponding to paragraph 6.a. 
could be obtained by assuming that the ply continues to carry failure 
loads (Path BA,) but can sustain no additional loads. This will corre- 
spond to an elastic perfectly plastic material. 

Ce Between the two extremes discussed in paragraphs 6.a. and 


6.b., other paths could result depending upon the nature of unloading of 


the affected ply. 


STRESS 


STRAIN 


Figure 1. Laminate Response Subsequent 
to Initial Failure 
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d. The affected ply may unload gradually. This would corre- 


spond to the path BA,- 


In the NOLAST program, the unloading option of the affected ply corresponds 
to the sath BB,A,- Introduction of other options into the program is not a 
problem. This is discussed in Section IV. 

The structure of the program 'NOLAST' is show in Figure 2. The main 
program calls Subroutines INPUT and SPLIN1 only once. Thereafter for every 


increment of loads, it calls Subroutines ELCON, ITER and OUTPUT in succession. 
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SECTION III 


USER'S INSTRUCTIONS FOR DATA PREPARATION 


The follewing sequence of punched cards will be necessary for 


conducting strength analysis of laminates: 


Card 1 HED CARD (12A6) 
Col 1-72 Any alphanumeric information 

| Card 2 MATERIALS (15) 

- | Col 1-5 MATYPE - number of material systems, (maximum = 5) 
Cards 3 NUMBER OF NODAL POINTS FOR STRESS-STRAIN PLOTS (715) 
Col 1-5 NEPOT - 0° tensile 


| 6-10 NEPOC - 0° compressive 
11-15 NEP90T 90° tensile 
| 16-20 NEP9OC 90° compressive 


! 21-25 NIP12 £- shear 
26-30 NPNUT - tensile Poisson's ratio 
| 31-35 NPNUC’ - Compressive Poisson's ratio 


Note 1. The above nodal points data includes the fictitious nodal 
point A (Figure 3). This is done to improve the accuracy 
of the curves in the last segment. Maximum number of 
nodal points for a stress-strain plot is limited to 30. 


| Cards 4 MATERIAL PROPERTY DATA (6F10.0) 
; Col 1-10 Strain Repeat until all stress-strain 
11-20 Stress data corresponding to 0° tension, 
21-30 Strain 0° compression, 90° tension, 90° 
31-40 Stress compression, shear, tensile Poisson's 
41-50 Strain ratio and compression Poisson's ratio 
51-60 Stress in sequence is accounted for. 
Cards 5 INITIAL ELASTIC PROPERTIES (6F10.0) 
Col 1-10 E1T Elastic modulus of 0° lamina in tension 
11-20 += £EIC Elastic modulus of 0° lamina in compression 
21-30 £E2T Elastic modulus of 90° lamina in tension 
31-40 E2C Elastic modulus of 90° lamina in compression 
41-50 G12 Shear modulis of 0° lamina 
Note 2. One set of (Cards 3 to 5) this data is required for each 
material system. 
H 


STRAIN 


Figure 3. Fictitious Nodal Point A 


Figure 4. Transformation of Axes 


ee ee 
pp ee a op Ra tn 


Card 6 
Col 


Cards 7 
Col 


Cards 8 
Col 


Cards 9 
Col 


Cards 10 
Col 


NUMBER OF LAMINATES (1415) 


1-5 NUMLAM - number of laminates 

CONTROL CARDS FOR EACH LAMINATE (1415) 

1-5 LAMINA - number of plies in a laminate 
(maximum = 15) 

6-10 NN - number of load combinations 
(maximum = 10) 

11-15 MLT -]| Exponents in Equations 5 to 7. 

16-20 MLC -| Until further experimentation, 

21-25 MIT ~| they may be assumed to be unity. 

26-30 MIC - 

31-35 MSH = 


36-40 NPRINT - Print output for each NPRINT increment (e.g. 
if NPRINT = 3, the output would be printed 
3rd, 6th, 9th etc. increments) except the 
output for first increment and every incre- 
ment after initial failure. 

Unloading option. It is to be taken 

equal to unity until additional options 

of unloading are added. For NOPSHN = 1, 
the affected ply or plies are unloaded 

in one step and loads transferred to the 
remaining plies. 


41-45 NOPSHN 


BOUNDARIES OF PLIES OF THE LAMINATE (6F10.0) 


1-10 = H(1) These are the distances of the boundaries 
11-20 = H(2) of the plies from a reference line as 
21-30 4H(3) defined in Figure 5. 

31-40 H(4) 


e 
H(LAMINA + 1) 


ORIENTATION OF PLIES (6F10.0) 


1-10 TH(1) Orientation angle in degrees of the 
first ply (Figure 4) 
11-20 TH(2) Orientation angle in degrees of the 


second ply (Figure 4) 
TH (LAMINA) Orientation angle in degrees of the 
last ply (Figure 4) 


MATERIAL TYPES ASSOCIATED WITH PLIES (1415) 


1-5 MAT(1) Material number of the first ply 
6-10 MAT (2) Material number of the second ply 
e 


e 
MAT(LAMINA) Material number of the last ply 
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Cards 11 


Col 


Note 


Note 


Note 


Note 


3.1 


3.2 


3.3 


4. 


LOAD CONDITIONS (3F10. n) 


1-10 A1(1) Stress resultant ANX in x-direction 
(Figure 4) for the first load condition. 

11-20 = A2(1) Stress resultant ANY in y-direction 
(Figure 4) for the first load condition. 

21-30 = A3(1) Stress resultant ANXY in xy-direction 


(Figure 4) for the first load condition. 
Repeat for all the number of load conditions NN of the Card 7. 


The components of stress resultants are computed by estimating 
strengths (NX, NY, and NXY) of the laminate and dividing the 
same by the number of desired increments. Experience indi- 
cates that 10-20 increments are adequate. 


ANX, ANY, ANXY of Cards 11 represent the maximum size of the 
increments. In the program the increment size get reduced a8 
the failure points based upon strains are approached. 


Prepare other sets of cards from Card 7 to Card 11 for other 
laminates defined by Card 6. 
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0.0 
Reference Line 


H(LAMINA + 1) 


Figure 5. Distance3 of Boundaries 
| of Plies of a Laminate 
from a Reference Line 
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SECTION IV 


PROGRAM MODIFICATIONS AND SAMPLE PROBLEMS \ 


le MODIFICATIONS 

NOLAST is the experimental computer program. It is being made 
available to designated recipients in the interest of timely exchange of 
technical information. Therefore user's experience communicated to 
AFFDL/FBR would be valuable for improving the program and would be 
appreciated. In the meantime, it is likely that the program would be 
modified by the users to suit their requirements. To help them in making 
modifications, the following comments are pertinent: 

a. The program, in the present form, allows for five material 
systems (isotropic and anisotropic), thirty nodal points (including the 
fictitious one) of stress-strain curves, fifteen plies in a laminate, and 
ten combinations of loads for each of the laminates. These maxima can be 
modified by changing the dimensions of the program. 

b. The program incorporates the criterion represented by 
Equations 5 to 7. Any other criterion, with constraints if any, may be 
introduced by insertion of suitable Fortran Statements in the Subrout ine 
OUTPUT (cards 118 to 250). 

c. Additional unloading options may be added to cards 373 and 
374 by assigning a suitable number to NOPSHN on card 297 in the Subroutine 
OUTPUT. 


d. For additional changes reference may be made to AFFDL/ FBR. 
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2. SAMPLE PROBLEMS 

In the present form, NOLAST is an in-core program requiring less 
than 60K octal. It was used to generate stress-strain responses of the 
following laminates reported herein: 

a.  (0/+45/90) | laminate with all the plies of one material and 
subjected to two combinations of loads. 

b. (0/99), . laminate with all the plies of one material subjected 
to one load combination. 

é. (0/0/+45/90) | laminate with (0/+45/90) plies of one material 


and 0° ply of the second material subjected to one load combination. 


Generation of the responses of the laminates (a, b, and c) required compila- 
tion and execution times of 8.28 seconds and 4.13 seconds respectively, on 


CYBER 7400 computer. 


The stress-strain data of both materials appears on pages 49 and 
50. The relevant inpuc information for the three laminates (a, b, and c) 
is printed on pages 51, 58, and 63 respectively. 

The output information on each of the printed pages consists of 
the following: 

(1) Load Increment, Load Combination, Laminate No. and 
NOPSHN No. 

(11) Average laminate elastic constants (EX, EY, UXY, GXY, 
ETAL and ETA2) during the load increment. 


(411) Components of the laminate strain increment. 
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(iv) Components of the laminate strain, stress resultant 
and stress. 

(v) Components of biaxial lamina strains and stresses. 

(vi) Components of equivalent lamina strains and corresponding 
stresses. 

(vii) Total energies due only to longitudinal tension and 
compression, transverse tenr.ion and compression, and shear for each of 
the material systems. This information is printed only for the first 
load increment. 

(viii) Energy ratios (longitudinal, transverse and shear), their 
gums for all plies and the sum of sums. 

(ix) Moduli of elasticity of all plies at the end of the 
load increment. 

(x) Information about failing plies and their modes. 

(xi) At the time of the failure of a ply or plies (e.g. 
page 53, Load Increment 12, Load Combination 1, Laminate No. 1, NOPSHN 
No. 1) additional information consists of the moduli of elasticity of 
all the plies and stress resultants of the affected plies during unloading 
and iterative process. The unloading process was indicated by the total 
contribution (1.010665) of the fourth ply exceeding unity. The present 
program does not have a provision to refine the laminate stresses corre~ 
sponding to the total contribution (1.010665) to those pertaining to the 
total contribution of unity. This can easily be done by interpolating 


the results of the current load increment and the previous one. These 


17 


remarks also apply to the final failure of the laminate. 

(xii) Final failure of the laminate (e.g. page 55, Load 
Increment 18, Load Combination 1, Laminate No. 1, NOPSHN No. 1) occurs 
when the stress resultants of the failing plies cannot be imposed upon 
the intact plies without being reduced or witsout the stiffness matrix 
A becoming singular. 


(xiii) Additfonal information about c, of the plies and A 


j ij 
of the laminate may be obtained by deleting the letter C in the beginning 


of the card numbers 245 and 247 of the Subroutine OUTPUT. 


APPENDIX A 


FORTRAN LISTING OF COMPUTER 


PROGRAM (NOLAST) 
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